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^ ■ Abstract 

O 
O , 

By encoding a qudit in a harmonic oscillator and investigating the d — > oo 
\ limit, we give an entirely new realization of continuous- variable quantum 

computation. The generalized Pauli group is generated by number and phase 
• operators for harmonic oscillators. 



The use of continuous-variable (CV) quantum computing allows information to be en- 



coded and processed much more compactly and efficiently than with discrete-variable (qubit) 
computing. With CV realizations, one can perform quantum information processes using 
fewer coupled quantum systems: a considerable advantage for the experimental realization 
^ ■ of quantum computing. The rapidly developing field of CV quantum information theory 
has applications to quantum error correction 0], quantum cryptography H and quantum 
teleportation ||, including an experimental realization of CV quantum teleportation [|J. 

At present, the proposed realization of CV quantum computation employs position eigen- 
states as a computational basis [|1]]; these states are approximated experimentally using 
highly squeezed states 0. Here, we introduce new CV realizations, where the generalized 
Pauli group is generated by the number operator N and a phase operator 8, and the compu- 
tational basis is given either by harmonic oscillator number eigenstates or phase eigenstates. 
These realizations are obtained formally by taking the d — > oo limit of the qudit, the d— 
dimensional generalization of the qubit, and are important for four key reasons: 

1. these CV realizations are entirely distinct from the position eigenstate computational 
basis realization, both in terms of the computational basis and in terms of the SUM 
gate; 

2. the SUM gate employs a standard Kerr optical nonlinearity to couple two modes; 
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3. these realizations give natural extensions of the qubit-based (discrete-variable) Pauli 
group, with a well-defined limiting procedure; and 

4. these realizations give a new implementation of the well-studied phase operator ||. 

In this letter, we review the generalized Pauli group for qudits, and we construct the 
generators of the generalized Pauli group in d dimensions using operators that will be shown 
to be expressible in terms of the SU(2) angular momentum and phase operators. This 
construction allows us to conveniently view the <i-dimensional space of qudits as the Hilbert 
space of a <i-dimensional irrep of SU(2). We also express the qudits in terms of harmonic 
oscillator states and investigate the d — > oo limit, and show that it is not the common 
generalization of the Pauli group for CV quantum information (i.e., the Heisenberg-Weyl 
group) with position eigenstates as the computational basis. Instead, we obtain a new CV 
realization, where the generalized Pauli group is generated by the number operator N and a 
phase operator 9. We also construct a second realization of qudits in terms of phase states 
and show that this realization is "dual" to the first realization given here. Finally, we discuss 
a realization of CV quantum computation in coupled harmonic oscillators. We establish a 
SUM gate, which serves as the CV analogue of the CNOT gate; this SUM gate employs a 
optical nonlinearity and is distinct in operation from the SUM gate suggested [1] for 
the position-eigenstate computational basis. 

We begin by reviewing the Pauli group of a qubit, and its generalization to the qudit. A 
qubit is realized as a state in a two-dimensional Hilbert space Ti.2- It is customary to choose 
two normalized orthogonal states, |0) and |1), to serve as a computational basis for 7i 2 . The 
unitary operators {X 2 = cr x , Z 2 = a z }, where the <7j are the Pauli spin matrices, generate 
the Pauli group using matrix multiplication. The elements of this group are known as Pauli 
operators and provide a basis of unitary operators on 7i 2 - 

Just as a qubit is realized as a state in a Hilbert space of dimension two, a qudit is realized 
as a state in a <i-dimensional Hilbert space Tid- It is useful to choose a computational basis 
{\s); s = 0,1,..., d— 1} for 7i dl which serves as the generalization of the binary basis 
{|0), |1)} of the qubit. 

A basis for unitary operators on Tid is given by the generalized Pauli operators ||[7|,[l| 

(X d ) a (Z d ) b , a,be0,l,...d-l, (1) 

where Xd and Z d are defined by their action on the computational basis as follows: 

X d \s) = \s + 1 (mod d)} , (2) 
Zd\s) = exp(27ris/<i)|s) . (3) 

The operators X d and Zd generate a group under matrix multiplication, known as the 
generalized Pauli group. Note that Xd and Z d are non-commutative and obey 

Z d X d = exp(2ni/d)X d Z d . (4) 

In the following, we give a representation of X d and Z d in the ci-dimensional Hilbert space of 
a SU(2) irrep of highest weight (angular momentum) j = (d — l)/2. The relevant generalized 
Pauli operators can be viewed in terms of SU(2) angular momentum and phase operators ||. 
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Consider the standard basis for the su(2) algebra { J z , j± = J x ± i j y }. Let {\j, m) z ] m = 
—j, . . . , j} denote the standard weight basis for the Hilbert space TL d =2j+i for an SU(2) irrep 
of highest weight (angular momentum) j. We use the simplifying notation of Vourdas || 
where we allow m to take all the integer (or half-integer) values modulo 2j ' + 1, thus defining 

\j,j + l)z = \j, ~j)z- 

With the computational basis defined to be 

\s) = \j,j - s) z , s = 0, 1, . . . ,d- 1 , (5) 

we now write the generators of the generalized Pauli group in terms of operators that act in 
a natural way on our SU(2) basis states. Because the basis states are eigenstates of J z , we 
have 

Z d i-» exp(2vri(j - J z )/d) , (6) 

which is unitary and satisfies Eq. (|3|). For the generalized Pauli operator X d , we use 

j 

X d i-> ^ \3i m )z(j,m + M- ( 7 ) 

m=-j 

One can easily check that X d satisfies Eq. (0) under the identification of Eq. (|7|) and that it 
is unitary. The operators X d and Z d satisfy Eq. (|) and together generate a representation 
of the generalized Pauli group for a qudit. It is convenient to view X d as the exponent of 
a Hermitian operator X d = exp(27ri 9 z /d), just as Z is generated by the operator J z . The 
operator 6 Z is the SU(2) phase operator of Vourdas 0. 

It is also possible to realize the operators X d and Z d as operators that act naturally on the 
space HLj of dimension d spanned by harmonic oscillator states of no more than d — 1 bosons. 
We define the computational basis to be the set of harmonic oscillator energy eigenstates 

\s) = \n = s)no, s = 0, 1, . . . ,d - 1 , (8) 

where N\n)n = n|n)no- Again, we apply the cyclic notation \d) = |0). Encoding a 
qudit in an oscillator is important not only for investigating the d —>■ oo limit, but also 
for realizing a qudit experimentally and for creating error-correcting codes for qudit-based 
computation [|1|]. 

In this Hilbert space, the generators X d and Z d of the generalized Pauli group for a qudit 
become 

d-l 

Xd^^2\s + l)(s\, Z d ^exp(2mN/d), (9) 

which are unitary on H. d . It is convenient to view X d as the exponent of a Hermitian operator 
9 Z , such that X d = exp(27ri z /d); the operator 9 Z is the Pegg-Barnett phase operator p. 
We will call this representation of the generalized Pauli group the number representation. 

The explicit realization of X d and Z d as unitary operators on the harmonic oscillator 
Hilbert space enables us to investigate the d — > oo limit; the limiting procedure for phase 
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operators has been thoroughly investigated ||,[|. In this limit, the computational basis 
remains the harmonic oscillator energy eigenstates (now including all states s = 0, 1, . . . , oo), 
following Eq. @. It is natural to generalize the operator Xd to a continuous transformation 
X(x), generated by the phase operator 9 Z ; i.e., 

X(x) = exp(ix§ z ), xeR. (10) 

Similarly, the CV generalization of Zd is obtained by replacing the finite angle 2ir/d in the 
expression for Zd in Eq. (J^j by the continuous angle z G R, so that we now have the unitary 
transformation Z(z) defined by 

Z(z) = exp(izN) . (11) 

Note that, in extending the above representation of the generalized Pauli group from qu- 
dits to CV representations, we do not obtain the usual generalization as the Heisenberg-Weyl 
group, with position x and momentum p operators as generators. Instead, the generalized 
Pauli operators are generated by the number operator N and the phase operator 6 Z ; these 
operators are in a sense "conjugate" like momentum and position, but there exist challeng- 
ing problems with defining the phase operator on the infinite-dimensional Hilbert space 
of the harmonic oscillator |5],|9|]. It is also interesting to note that the states of the computa- 
tional basis for the limiting case remain harmonic oscillator energy eigenstates, not position 
(or momentum) eigenstates or squeezed Gaussians as are commonly used for CV quantum 
computing. In what follows, we will refer to this representation as the number representation 
of the generalized CV Pauli group. 

It is possible to construct another realization of Xd and Zd in the Hilbert space Tid for an 
irrep of SU(2) where the computational basis is given by SU(2) phase states. This represen- 
tation is "dual" to the number representation. Consider the relation iX 2 = exp(i(7r/2)X2) 
for qubits; i.e., that 

|1) =X 2 |0) = (-i)e i( " /2)X2 |0). (12) 

The Pauli operator X 2 has two interpretations, each of which can be generalized in a different 
way. In the number representation, we interpreted X2 as a number state raising operator 
|1) = X 2 |0) and generalized this operator as such. However, using the relation (|12D , we can 
also view X 2 as a rotation. (Using the su(2) representation X 2 = 2J X , this rotation is about 
the x-axis.) Thus, the state |1) is obtained (up to a phase) by rotating |0) by an angle tt 
about the x-axis. The computational basis states needed for this type of generalization to 
qudits are "SU(2) phase states" and have been investigated by Vourdas || (although using 
rotations generated by J z rather than J x ). These states form an orthonormal basis for the 
SU(2) irrep and are "dual" to the usual weight basis. 

Let {\j, m) x ; m = —j, . . . , j} be the weight basis for an SU(2) irrep of angular momentum 
j = (d — l)/2, where J x rather than J z is diagonal; i.e., J x \j, m) x = m\j,m) x . For this 
representation, we define the computational basis states to be 

1 

\s) = —= \ exp(2irims/d)\j 1 m) x , d odd, (13) 
yd . 

m=-j 

1 j 1 

= ex P( 27ri ( m + 2) s / d )\h m )x , deven. (14) 
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These states form an orthonormal basis for Tid @. They are referred to as SU(2) phase 
states because they are eigenstates of the SU(2) phase operator, defined below. 
The generalized Pauli operator Xd on this computational basis is given by 

Xd i— > exp(27ri J x /d) , d odd, (15) 
X d i— > exp(— m/d) exp(27riJ x /d) , eleven, (16) 

satisfying Eq. (0). Note that (Xd) d = 1 for both j integral and half-integral (i.e., spinor). 
The generalized Pauli operator Zd is given by 

d-l 

Z d ^Ys e M^is/d)\s)(s\, (17) 

s=0 

which is unitary and satisfies Eq. (0). Note that we can express Zd as the exponent of a 
Hermitian operator, 

d-l 

Z d = exp{2m6 x /d), § x = ^ s\s) (s\ ; (18) 

s=0 

the operator Q x is the SU(2) phase operator. 

Note that this representation of the generalized Pauli group is "dual" to the number 
representation in the same sense that the position and momentum representations of the 
harmonic oscillator are dual. For the number representation, the computational basis states 
are eigenstates of J z , and the phase operator 9 Z generates the "ladder" transformations. 
In the phase representation given here, the computational basis states are eigenstates of 
the phase operator 9 X , i.e., "phase eigenstates", and it is J x which generates the ladder 
transformations via rotations about the x-axis. Both of these representations can be con- 
sidered natural generalizations of the qubit case, because the standard computational basis 
|0) = | \) z and |1) = ||, —\)z are both eigenstates of J z and phase eigenstates of 6 X . 

As with the number representation, this phase representation of the generalized Pauli 
group can be expressed in a harmonic oscillator Hilbert space. Again, the d — > oo limit 
yields challenging problems: it is well known that phase eigenstates do not exist in the 
infinite-dimensional Hilbert space of the harmonic oscillator ||. 

In any experimental realization, the problems associated with taking the d —>■ oo limit 
would not arise. A physically realistic system would have a finite energy cutoff (and an 
associated resolution in time and thus phase), and so experimental CV computation would 
in actuality involve qudits with finite (although possibly very large) d. As a result of our 
well-defined limiting procedure for qudits, the above realization of CV quantum computation 
is applicable to such a physically realistic system. 



TO perform universal CV computation JT0[, it is necessary to be able to realize an 



arbitrary unitary transformation on a single qudit, and to have a controlled two-qudit in- 
teraction gate such as the SUM gate Considering an optical realization, an arbitrary 
unitary transformation on a single qudit, to any desired precision, can be performed effi- 
ciently using a combination of linear optics, parametric down-conversion, and a nonlinear 
optical Kerr medium pi . By this combination, one can approximate (to arbitrary accuracy) 
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any polynomial Hamiltonian in fit and a. Of particular importance is to realize the Fourier 
transform operation on a single qudit, which takes number eigenstates to phase eigenstates 
and vice versa. This operation is the generalization of the Hadamard transformation for 
qubits. 

For quantum computation, we must also realize a gate that performs a two-qudit inter- 
action. Consider two oscillators coupled by the four-wave mixing interaction Hamiltonian 
XN1N2 = x^-i^i^l^- This Hamiltonian for an optical system describes a four- wave mix- 
ing process in which \ is proportional to the third-order nonlinear susceptibility [[□]]. Let 



oscillator 1 be in a state |sx)i encoded in the number state basis, and let oscillator 2 be in 
a state ^2)2 encoded in the phase state basis. This interaction Hamiltonian generates the 
transformation 

e -m*\ Sl ) x ® | S2 ) 2 = ® \A Sl + S2 ) 2 . (19) 

Z7T 

Thus, with time t = 2ttx~ 1 , this Hamiltonian generates the SUM transformation |si)i (g) 

In summary, we have presented a new form of continuous variable computation in terms of 
number and phase operators. This new approach has the advantage over position-eigenstate 
CV computation in that the computational basis states, for large but finite d, are well- 
defined and obtainable, and do not require "infinite-squeezing" of Gaussian wavepackets. 

This project has been supported by an Australian Research Council Large Grant and by 
a Macquarie University Research Grant. We acknowledge helpful discussions with Samuel 
Braunstein. 
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